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Digital Cancelation of Self-Interference for Single-Frequency Full-Duplex 
Relay Stations via Sampled-Data Control 


Hampei Sasahara*, Masaaki Nagahara*, Kazunori Hayashi*, and Yutaka Yamamoto* 

Abstract: In this article, we propose sampled-data design of digital filters that cancel the continuous-time effect of 
coupling waves in a single-frequency full-duplex relay station. In this study, we model a relay station as a continuous¬ 
time system while conventional researches treat it as a discrete-time system. For a continuous-time model, we propose 
digital feedback canceler based on the sampled-data H°° control theory to cancel coupling waves taking intersample 
behavior into account. We also propose robust control against unknown multipath interference. Simulation results are 
shown to illustrate the effectiveness of the proposed method. 
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1. Introduction 

In wireless communications, relay stations have been used 
to relay radio signals between radio stations that cannot di¬ 
rectly communicate with each other due to the signal attenu¬ 
ation. More recently, relay stations are used to achieve a cer¬ 
tain spatial diversity called cooperative diversity to cope with 
fading channels [1]. On the other hand, it is important to ef¬ 
ficiently utilize the scarce bandwidth due to the limitation of 
frequency resources [2], while conventional relay stations com¬ 
monly use different wireless resources, such as frequency, time 
and code, for their reception and transmission of the signals. 
For this reason, a single-frequency full-duplex relay station, in 
which signals with the same carrier frequency are received and 
transmitted simultaneously, is considered as one of key tech¬ 
nologies in the fifth generation (5G) mobile communications 
systems [3]. In order to realize such full-duplex relay stations, 
self-interference caused by coupling waves is the key issue [4], 

Fig. 1 illustrates self-interference by coupling waves. In this 
figure, radio signals with carrier frequency / are transmitted 
from the base station (denoted by BS). One terminal (denoted 
by Tl) directly receives the signal from the base station, but the 
other terminal (denoted by T2) is so far from the base station 
that they cannot communicate directly. Therefore, a relay sta¬ 
tion (denoted by RS) is attached between them to relay radio 
signals. Then, radio signals with the same carrier frequency / 
are transmitted from RS to T2, but also they are fed back to 
the receiving antenna directly or through reflection objects. As 
a result, self-interference is caused in the relay station, which 
may deteriorate the quality of communication and, even worse, 
may destabilize the closed-loop system. 

To tackle with the issue of self-interference cancelation, 
many methods have been proposed for single-frequency full- 
duplex systems. Analog cancelation has been proposed in 
[5],[6], in which analog devices are used for canceling cou¬ 
pling waves. Since coupling wave paths are physically ana¬ 
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Fig. 1 Self-interference 


log systems and there is no quantization problem, this design is 
theoretically the most ideal except for implementation issues. 
On the other hand, digital cancelation has attracted increasing 
attention, in which the interference is subtracted in the digital 
domain by using digital signal processing techniques [7]—[13]. 
Digital cancelers benefit easy implementation on digital devices 
in exchange for the response between sampling instants. In ad¬ 
dition, spatial domain techniques, called antenna cancelation, 
has been also proposed in [4],[14], in which they try to reduce 
the interference by arranging antenna placement. See [4],[8] 
for details. 

For the problem of self-interference, a pre-nulling method 
L10] and adaptive methods [7],[12] have been proposed to can¬ 
cel the effect of coupling waves. In these studies, a relay station 
is modeled by a discrete-time system, and the performance is 
optimized in the discrete-time domain. However, radio waves 
are in nature continuous-time signals and hence the perfor¬ 
mance should be discussed in the continuous-time domain. In 
other words, one should take account of intersample behavior 
for coupling wave cancelation. 

In theory, if the signals are completely band-limited below 
the Nyquist frequency, then the intersample behavior can be re¬ 
stored from the sampled-data in principle [15], and the discrete¬ 
time domain approaches might work well. However, the as¬ 
sumption of perfect band limitedness is hardly satisfied in real 
signals since 
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1. real baseband signals are not fully band-limited, 

2. pulse-shaping filters, such as raised-cosine filters, do not 
act perfectly, 

3. and the nonlinearity in electric circuits adds frequency 
components beyond the Nyquist frequency. 

One might think that if the sampling frequency is fast enough, 
then the assumption is almost satisfied and there is no problem. 
But this is not true; firstly, the sampling frequency cannot be 
arbitrarily increased in real systems, and secondly, even though 
the sampling is quite fast, intersample oscillations may happen 
in feedback systems [16, Sect. 7]. 

To solve the problem mentioned above, we propose a new 
design method for coupling wave cancelation based on the 
samplecl-data control theory [16],[17]. We model the transmit¬ 
ted radio signals and coupling waves as continuous-time sig¬ 
nals, and optimize the worst case continuous-time error due to 
coupling waves by a digital canceler. This is formulated as a 
sampled-data H°° optimal control problem, which can be solved 
via the fast-sampling fast-hold (FSFH) method [18],[19]. We 
also propose robust feedback cancelers that can take account of 
uncertainties in coupling wave path characteristic such as un¬ 
known multipath interference due to, for example, large struc¬ 
tures that reflect radio waves, or the change of weather con¬ 
ditions [20]. Design examples are shown to illustrate the pro¬ 
posed methods. 

The present manuscript expands on our recent conference 
contributions [21],[22] by incorporating robust feedback con¬ 
trol into the formulation. 

The remainder of this article is organized as follows. In Sec¬ 
tion 2, we derive a mathematical model of a relay station con¬ 
sidered in this study. In Section 3, we propose sampled-data 
H°° control for cancelation of self-interference. Here we also 
discuss robust control against uncertainty in the delay time. In 
Section 4, simulation results are shown to illustrate the effec¬ 
tiveness of the proposed method. In Section 5, we offer con¬ 
cluding remarks. 

Notation 

Throughout this article, we use the following notation. We 
denote by I? the Lebesgue space consisting of all square inte- 
grable real functions on [0, oo) endowed with L 2 norm || • H 2 - 
The symbol t denotes the argument of time, .y the argument of 
Laplace transform and z the argument of Z transform. These 
symbols are used to indicate whether a signal or a system is of 
continuous-time or discrete-time. The operator e~ Ls with non¬ 
negative real number L denotes the continuous-time delay op¬ 
erator with delay time L. For a matrix A, cr(A) denotes the 
maximum singular value of A. 

2. Relay Station Model 

In this section, we provide a mathematical model of a relay 
station with self-interference phenomenon. 

Fig. 2 depicts a single-frequency full-duplex relay station im¬ 
plemented with a digital canceler [6]. A radio wave with carrier 
frequency / from a base station is accepted at the receiving an¬ 
tenna and amplified by the low noise amplifier (LNA). Then, 
the received signal is demodulated to a baseband signal by the 
demodulator, and converted to a digital signal by the analog- 
to-digital converter (ADC). The obtained digital signal is then 
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Fig. 3 Simple Block Diagram of Relay Station 


processed by the digital signal processor (DSP) into another 
digital signal, which is converted to an analog signal by the 
digital-to-analog converter (DAC). Finally, the analog signal is 
modulated to a radio wave with carrier frequency /, amplified 
by the power amplifier (PA) and transmitted by the transmis¬ 
sion antenna. A problem here is that the transmitted signal will 
again reach the receiving antenna. This is called coupling wave 
and causes self-interference, which deteriorates the communi¬ 
cation quality. 

Fig. 3 shows a simplified block diagram of the relay station. 
In Fig. 2, we model LNA and PA in Fig. 2 as static gains, a\ 
and a 2 , respectively. The modulator is denoted by A! and the 
demodulator by D. We assume that the coupling wave channel 
is a flat fading channel, that is, all frequency components of a 
signal through this channel experience the same magnitude fad¬ 
ing. Then the channel can be treated as an all-pass system. In 
this study, we adopt a delay system, re ~ Ls , as a channel model, 
where r > 0 is the attenuation rate and L > 0 is a delay time. 
The block named “Digital System” includes ADC, DSP and 
DAC in Fig. 2. 

In this article, we consider the quadrature amplitude modu¬ 
lation (QAM), which is used widely in digital communication 
systems, as a modulation method. QAM transforms a trans¬ 
mission signal into two orthogonal carrier waves, that is a sine 
wave and a cosine wave. We assume the transmission signal 
u(t) is given by 
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( 1 ) 


In this expression, g(t) is a general pulse-shaping function, li is 
the sampling period, and u[, u@ denote respectively the in-phase 
and the quadrature components of a transmission symbol. We 
assume that the support of the Fourier transform G(ju>) of g(t) is 
finite and the bandwidth is much less than 4 nf. In other words, 
there exists a frequency f g (0 < f g <£ /) such that \G(ja>)\ = 0 
for any a> i (-2nf g , 2itf g ). Then the modulated signal u(t) can 
be written as [23, Chap. 2] 



Fig. 4 Structure of a demodulator 


u(t)=Mu(t) 

= ^ g(t - kh)(u[ cos 2 trft - u® sin 2nft). (2) 

k 

On the other hand, the demodulation operator D is a linear 
operator satisfying DM = 1 [23]. Fig. 4 shows the block dia¬ 
gram of D. In this block diagram, H [d (jaj) is the ideal low-pass 
filter with cut-off frequency f c satisfying 


HAM = 


1, if a) < 2 7if c , 
0, otherwise. 


(3) 


The cut-off frequency is chosen to satisfy f g <sc f c /. By the 
linearity of D , we obtain an equivalent block diagram shown in 
Fig. 5 to Fig. 3. Here 


u(t - L)=^ g(t - L - k/i)|(a[ cos(2 nfL) 
k 
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Fig. 5 Equivalent Block Diagram of Relay Station 



Fig. 6 Relay Station Model 


Fig. 7 Digital System 



+(n[ sin(2;r/L) - u g cos(2 nfL)) sin(2;r/f)j. (4) 
Thus, we have 

cos(2 nft) ■ u(t - L) 

= \y\g(t- L - kh)l u[ cos(2 nfL) + uf sin(2 nfL) 

Z k 1 

+(u[ cos(2 nfL) + uf sin(2 nfL)) cos(4;r/f) 

+(u[ sin(2 nfL) - u® cos(2 nfL)) sin(4^/?)}. (5) 

From this, we have 


Finally we obtain a relay station model depicted in Fig. 6. 
By this figure, we can see that the relay station with self¬ 
interference is a feedback system. In practice, the gain of PA in 
Fig. 2 is very high (e.g., a-i = 1000) and the loop gain becomes 
much larger than 1, and hence we should discuss the stability as 
well as self-interference cancelation. To achieve these require¬ 
ments, we design the digital controller in the digital system, 
which is precisely shown in Fig. 7. 

In Fig. 7, ADC in Fig. 2 is modeled by an anti-aliasing analog 
filter F(s) with an ideal sampler Si, with sampling period li > 0, 
defined by 


2HAcos(2ttft) ■ u(t - L)] 

= g(t - L - kh){u[ cos(2 nfL) + u k sin(2;r/L)}. (6) 

k 

In the same way, we have 

2/7id[- sin(2nft) ■ u(t - L)] 

= ^ j g(t-L- kh){-u‘ k sin(2 nfL) + uf cos(2;r/L)}.(7) 

k 

Finally, we have the following relation: 

U]{t)=D(axrafuit - L)) 

=aA L u(t - L), (8) 


where a := a^r and 

cos(2 ttfL) sin(2 nfL) 

L ' - sin(2 nfL) cos(2 ttfL) 


S h : {>’o(01 •-» [ydln]} ■■ yd[n\ = yo (nh), 
n = 0 , 1 , 2 ,.... 

For the DSP block in Fig. 2, we assume a digital filter de¬ 
noted by K(z), which we design for self-interference cancela¬ 
tion. DAC in Fig. 2 is modeled by a zero-order hold, 7Y/,, de¬ 
fined by 

'H h : {u d [n ]} ^ [MoWJ : «o(0 = u d [n], 

t€[nh,(n+l)h),n = 0 , 1 , 2 ,..., ^ ' 

and a post analog low-pass filter denoted by P(s). We assume 
that F(s) and P(s) are proper, stable and real-rational transfer 
function matrices. Note that a strictly proper function is nor¬ 
mally used for F(s) and it is included in the assumption. 

3. Feedback Control 

Fig. 8 shows the block diagram of the feedback control sys¬ 
tem of the relay station. For this system, we find the digital 
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Fig. 8 Feedback Canceler 


controller, K(z), that stabilizes the feedback system and also 
minimize the effect of self-interference, z := v — u, for any v. To 
obtain a reasonable solution, we restrict the input continuous¬ 
time signal v to the following set 

WLr := {v = Ww : w € L 2 , \\w\\ 2 = 1}, (12) 

where W is a continuous-time LTI system with real-rational, 
stable, and strictly proper transfer function W(s). Under this 
assumption, we first solve a nominal control problem where 
all system parameters are previously known. Then we propose 
a robust controller design against uncertainty in the coupling 
wave paths. 

3.1 Nominal Controller Design 

Here we consider the nominal controller design problem for¬ 
mulated as follows: 

Problem 1 Find the digital controller (canceler) K(z) that sta¬ 
bilizes the feedback system in Fig. 8 and uniformly minimizes 
the L 2 norm of the error z = v - u for any v e WLr. 

This problem is reducible to a standard sampled-data H°° 
control problem [16],[17]. To see this, let us consider the block 
diagram shown in Fig. 9. Let T zw be the system from w to z. 
Then we have 


z = v - u = T zw w (13) 

and hence uniformly minimizing || z [| 2 for any v £ WLr is equiv¬ 
alent to minimizing the H°° norm of T zw , 

II77JL = sup ||r w w|| 2 . (14) 

weL 1 , |Mb = l 


Let 2(i) be a generalized plant given by 

L(.s) = 


W(s) - P(s ) 

F(s)W(j) ae~ Ls A L F(s)P(s) 


By using this, we have 

T zw (s) = T(Z(s)/H h K(z)S h ), 


(15) 


(16) 


where F denotes the linear-fractional transformation (LFT) 

[17] . Fig. 10 shows the block diagram of this LFT. Then 
our problem is to find a digital controller K(z) that minimizes 
HTT-mIIoo. This is a standard sampled-data H°° control prob¬ 
lem, and can be efficiently solved via FSFH approximation 

[18] ,[19],[24], 

Note that if there exists a controller K(z ) that minimizes 
HT^IU, then the feedback system is stable and the effect of 
self-interference z = v - u is bounded by the H°° norm. We 
summarize this as a proposition. 


Proposition 1 Assume ||7’ ; ., v || 00 < y with y > 0. Then the feed¬ 
back system shown in Fig. 8 is stable, and for any v £ WLr we 
have ||v - u \| 2 < y. 



Fig. 9 Block Diagram for Feedback Canceler Design 



Fig. 10 LFT T zw = F(F < H h KS h ) 


Proof First, if the feedback system is unstable, then the H°° 
norm becomes unbounded. Next, for v £ WLr there exists w £ 
Lr such that v = Ww and ||w|| 2 = 1. Then, inequality HT^Hoo < 
y gives 

llv - u\\ 2 = ||r w w|| 2 < HT-JUMk < y. (17) 


□ 

3.2 Robust Controller Design against Multipath Interfer¬ 
ence 

In practice, the characteristic of the coupling wave channel 
changes due to, for example, large structures that reflect radio 
waves. In this situation, it is difficult to predict the coupling 
wave paths beforehand, and hence there must be uncertainties 
in the paths. Under this uncertainty, the nominal controller 
may lead to deterioration of cancelation performance, and even 
worse, it may make the feedback system unstable. To solve 
this problem, we propose robust controller design against the 
uncertainty. 

Let us assume the characteristic of the coupling wave paths 
in Fig. 6 is perturbed as 

M 

re~ Ls i — > re~ Ls + r ie ~ LiS , (18) 

i=i 

where r, and L, are the attenuation ratio and the delay time of 
the r'-th path, respectively. Note that M represents the number 
of additional paths. Since the additional paths are detour paths, 
we assume 

Li>L, i = 1,2,... ,M. (19) 

Then the characteristic of the feedback path in Fig. 8 is per¬ 
turbed as 

M 

ae~ Ls A L i — > ae~ Ls A L + ^ aie~ LiS Ai i (20) 

i=i 

where a,- := a ] U 2 r ! . Define the error transfer function matrix 
E(s) as 
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Fig. 11 Relay Station Model with Perturbation 
M 


E(s) := V -i e-^~ L)s A Lj . L . 

a 

i-i 


( 21 ) 


Since A L is a rotation matrix on R 2 and the angle is 2 nfL clock¬ 
wise, we have 


M 

ae~ Ls A L + ^ a t e~ LiS A Li 

i= 1 


M 

= ae~ Ls A L [l + £ ^e- (Li - L)s Al l A L ^j 
1=1 

= ae~ Ls A L (/ + E(s )). (22) 

Take a frequency weighting function matrix W 2 {s) that is real 
rational and satisfies 


Feedback canceller 



Fig. 12 Feedback Cancelation: input signal (dash-dot line), reconstructed 
signal it by feedback canceler (solid line) 

To solve this problem, we adopt the finite dimensional Q- 
parametrization where we limit feasible controllers [26]. Then, 
the constraints are represented by linear matrix inequalities 
(LMI’s) and the problem can be efficiently solved via numeri¬ 
cal optimization software such as SDPT3 or SeDuMi on MATLAB 
[27],[28]. For more details, see [26]. 


o-{E{ju>)) < cr(W 2 (jaj)), 


(23) 


for all weR. Since A /_/, is an orthogonal matrix, the equation 
(21) gives 


<r(E(M) < £ < £ £ 


(24) 


Then the uncertainty in the coupling wave paths can be mod¬ 
eled as multiplicative perturbation, that is, for any M > 0, 
Tj > 0, Lf > L (i = 1,..., M ), we have 


4. Design Examples 

In this section, we show simulation results to illustrate the 
effectiveness of the proposed methods. 

We assume that the sampling period h is normalized to 1, 
the carrier frequency / is 10000 Hz, the attenuation rate of the 
coupling wave channel r = 0.2, and the time delay L = 1. Note 
that sampling frequency is 1 Hz, which is much smaller than 
the carrier frequency. Note also that the time delay is equal to 
the sampling period h. We assume the low noise amplifier a\ = 
1. An anti-alias analog filter is not employed in this examples, 
namely we assume F(s) = I. The post filter P(s) is modeled by 


M 

ae~ Ls A L + ^ aie~ LiS A Li 
1=1 


e {ae~ Ls A L { 1 + A(s)W 2 {s )): ||A|L < 1}. 


(25) 


From the inequality (24), we can take 



(26) 


where e is an appropriately small and positive number. 

Based on the formulation of the uncertainty, we consider the 
block diagram shown in Fig. 11, where W\ (s) plays the same 
role as W(s) in the nominal controller design. Let T ZlWl be the 
system from w\ to zi and T Z1W1 be the system from W2 to z 2 - If 
ll^wilU is finite and ||r aW2 |U < 1. then the feedback system is 
robustly stable, that is, the feedback system is internally stable 
for all A satisfying HAH^ < 1 from the small gain theorem for 
sampled-data control systems [25]. 

Now we formulate the robust controller design problem as 
follows: 


Problem 2 Find the digital controller (canceler) K(z) that min¬ 
imizes ||7'- lV i. 1 ||oo subject to ||77 2H , 2 ||oo < 1. 


P(s) 


1 

0.001 y + 1 


I. 


We also assume the transmission gain to be 


a 2 = 1000, 


(27) 


(28) 


that is, 60 dB. The frequency characteristic W(s) is modeled by 

W(s) = —^—1. (29) 

2s + 1 

With these parameters, we compute the //"-optimal nominal 
controller K(z) by FSFH with discretization number N = 16. 

With this controller, we simulate coupling wave cancelation 
with a random rectangular wave input with period 4 s filtered by 
the low-pass filter P(s). Note that this signal contains frequency 
components beyond the Nyquist frequency, n //t = n [rad/sec], 
although the frequency of the wave, 7t/8/z = 7t/8 [rad/sec] is 
much lower than n. Fig. 12 shows the reconstructed signal u 
in the feedback system (see Fig. 8). The feedback system is 
guaranteed to be stable and the canceler achieves small recon¬ 
struction errors as shown in Fig. 13. 

Next, let us consider uncertainty in the coupling wave paths. 
If the characteristic of the coupling wave paths changes, then 
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Error 



Fig. 13 Coupling wave effect \v(t) - u(t)\ by feedback canceler shown in 
Fig. 8 


Robust canceller 



Fig. 15 Robust Cancelation: input signal (dash-dot line) and recon¬ 
structed signal (solid line) with the perturbation. 


Feedback canceller 



Fig. 14 Feedback Cancelation: input signal (dash-dot line) and recon¬ 
structed signal (solid line) with the perturbation 

the nominal controller may not work well. To see this, let us 
design the nominal controller with 

a 2 = 100, (30) 

that is, 40 dB, and the other parameters are the same as above. 
Then perturbing the nominal coupling wave path to be 

re~ Ls + r x e~ L ' s (31) 

where r\ = 0.07 r,L\ = 1.1L. It results in instability as shown 
in Fig. 14. 

To overcome this, we use the robust controller proposed in 
subsection 3.2, M = \, r\ = O.lr with VFi(s) = VT(^) given in 
(29) and W 2 (s) as in (26). The dimension of (9-parametrization 
is 8 with the FSFH number N = 4. Fig. 15 shows the recon¬ 
structed signal, by which we can observe the robust controller 
works well. 

5. Conclusions 

In this paper, we have proposed feedback controller design 
for self-interference cancelation in single-frequency full-duplex 
relay stations based on the sampled-data H°° control theory. In 
particular, we proposed robust controller design against the un¬ 
known additive multipath. Simulation results have been shown 


to illustrates the effectiveness of the proposed cancelers in view 
of stability and robust stability. Future work may include FIR 
(Finite Impulse Response) filter design and adaptive FIR filter¬ 
ing as discussed in [29],[30]. 
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